Temperature dependent spatial oscillations in the correlations of the XXZ spin chain 
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We study the correlation < a^a^ > for the XXZ chain 
in the massless attractive (ferromagnetic) region at positive 
temperatures by means of a numerical study of the quantum 
transfer matrix. We find that there is a range of temperature 
where the behavior of the correlation for large separations is 
oscillatory with an incommensurate period which depends on 
temperature. 
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I. INTRODUCTION 

Recently we Q- [|J have studied the finite temperature 
spin correlation function 

S z {n; T, A) = Tra z a z n e- H / kT /Tre- H / kT (1) 

for the spin 1/2 XXZ spin chain §]- |] 
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(2) 



In [Q we made a finite size study of S z (n; T, A) for chains 
of size up to L = 18 and found that for —1 < A < 
there is a range of temperature where the correlation 
function changes sign as a function of n starting from 
negative values for n — 1. For T sufficiently low the 
correlations for —1 < A < are negative for all n whereas 
for sufficiently large T the correlations are positive for 
all n. We referred to this change in sign as a quantum 
classical crossover. 

In both [jlj and Q we discussed the representation of 
this correlation function for the chain with L — > oo in 
terms of an expansion in terms of eigenvalues and eigen- 
vectors of the quantum transfer matrix Q 



S'(n;T,£) = J2M*i/>*) 



(3) 



where Ao is the maximum eigenvalue of the quantum 
transfer matrix and is directly related to the free energy 
per lattice site by 



/= -/cTlnAo + A/2 



(4) 



This quantum transfer matrix is defined as the limit 
N — > oo of a matrix of dimension 2 N where N, which 
is even, is the "Trotter number" (and is not to be con- 
fused with the chain length L which is taken to infinity.) 
We have numerically studied this expansion for values of 
N as large as 16 and have found several features of the 
correlations which have not previously been seen which 
substantially expands our understanding of the phenom- 
ena of quantum-classical cross over. 

From our study we find that as the temperature in- 
creases from zero the following phenomena occur: 

1. There is a temperature Ty(A) below which all am- 
plitudes Aj are negative and at which temperature 
many Aj vanish. There is no other temperature at 
which any Aj vanishes. For T < Ty(A) all eigen- 
values are real. This temperature vanishes as X/N 
as iV — > oo. We will refer to TV (A) as the universal 
vanishing temperature. 

2. There is a temperature Tc(A) > TV (A) above 
which real eigenvalues collide and become complex 
conjugate pairs. At the temperatures where the 
collision of the eigenvalues occurs the amplitudes 
Aj diverge as a square root. In the limit N —> oo 
we have T C (A) -> 2V (A). 

3. There is a temperature Tl(A) > Tc(A) where the 
eigenvalue in (^|) with the largest magnitude be- 
comes complex. As N — > oo T^(A) approaches 
a nonzero limiting value from above. We refer to 
Tj,(A) as the lower crossover temperature. 

4. There is a temperature TV (A) > TV (A) where the 
eigenvalue in (|3|) with the largest magnitude be- 
comes real and the corresponding amplitude is pos- 
itive. At TV (A) eigenvalues with different quantum 
numbers cross. As N — > oo TV (A) approaches its 
limiting value from below. We refer to TV (A) as 
the upper crossover temperature. 

We conclude that for T between T L (A) and TV (A) 
the correlation S z (n; T, A) oscillates and changes sign an 
infinite number of times as n — » oo. 

In the remainder of this letter we will present the evi- 
dence to support and illustrate these conclusions. 
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II. THE UNIVERSAL VANISHING 
TEMPERATURE TV (A) 

At T = it has been known for some time 
when n — > oo that for —1 < A < 



S"(n;0,A) 



1 



2 6>n 2 



+ (-1) 



(5) 



where cos7r6> = —A and the constant C(A) has recently 
been conjectured JTlJ] . In the very low temperature limit 
T — » and n — ► oo with nT = r fixed where conformal 
field theory is applicable [fL2| - [0] the large n behavior of 
S z (n,T, A) is obtained from (f3|) by the replacement 
n — » (nT/2)~ 1 sinIiKr/2 where n — tt(\ — 9)/ sin-Trf?. 

From our numerical study we find that for sufficiently 
low temperatures all eigenvalues Xj are real (although not 
all are positive) and all amplitudes Aj are negative. We 
illustrate this in table 1 where we show the first 13 eigen- 
values and amplitudes in (|J) for A = — .5,T = .1,N — 
12. We choose A = — .5 only for purposes of illustration. 
All phenomena discussed in this paper occur generally 
for —1 < A < 0. Here, and in the subsequent tables, we 
give only the eigenvalues in the subspace which is odd 
under spin inversion because the matrix elements in the 
complementary subspace are identically zero. We also 
give the quantum number k which plays for the quan- 
tum transfer matrix the role which momentum plays for 
the row transfer matrix. We note for N = 0(mod4) that 
k = 0, ±1, • • • , ±(N/4 - 1), N/4 and for N = 2(mod4) 
that 0, ±1, • ■ • , ±( A^ 2 -). The eigenvalues for ±k are de- 
generate because of reflection symmetry of the quantum 
transfer matrix. 



the first 13 eigenvalues and amplitudes for this case. Be- 
cause the temperature is slightly above the universal van- 
ishing temperature the sign of the amplitudes with the 
vanishingly small magnitudes are positive. The order of 
the eigenvalues 6-11 in table 2 is correct even though to 
the number of significant figures given the eigenvalues are 
identical. It should be noted that in table 2 all eigenval- 
ues are still real. 



j 


h 


IVAo| 


Aj phase 


\Aj\ 


Aj phase 


k 


1,2 


1,2 


.4791320 


0.0 


4.285003xl0~ 2 


7T 


±1 


3 


3 


.3177245 


7T 


7.130329xl0~ 3 


7T 





4,5 


4,5 


.2424384 


0.0 


5.180946xl0" 2 


7T 


±2 


6,7 


7,8 


.1737014 


0.0 


1.726434xl0~ 10 


0.0 


±2 


8,9 


9,10 


.1737014 


0.0 


5.754777xl0~ 1:L 


0.0 


±1 


10 


13 


.1737014 


0.0 


4.316086xl0~ n 


0.0 





11 


6 


.1737014 


0.0 


5.186948xl0" 2 


7T 


3 


12,13 


9,10 


.1328673 


7T 


4.055662xl0" 3 


7T 


±1 



TABLE II. Eigenvalues at Tv(A). The first 13 eigenvalues 
Aj and amplitudes Aj for A = -.5,T = 0.137832224 and 
N = 12. The value of A is 10.7928. The column ji is the 
value of j of the corresponding eigenvalue for T = .1 in table 
1. 

The numerically determined values of Ty{— .5) for N = 
10, 12, 14, 16 are given in table 6. The N dependence is 
extremely well fit for 6 < N < 16 by T V (-.5;N) = 
1 .6539 /A^ and therefore TV (A) vanishes as N — > oo. We 
thus conclude that the data for T < TV (A) and N finite 
is not relevant to the genuine quantum system with TV — > 
oo. 



j 


lA./Aol 


Aj phase 


14* | 


Aj phase 


k 


1,2 


.5617564 


0.0 


2.906227x10^ 


7T 


±1 


3 


.3954749 


7T 


5.619120xl0~ 3 


7T 





4,5 


.3346469 


0.0 


3.258793xl0" 2 


7T 


±2 


6 


.2323767 


0.0 


4.927463 xlO" 2 


7T 


3 


7,8 


.1911818 


0.0 


1.733959xl0~ 2 


7T 


±2 


9,10 


.1833028 


7T 


4.491119xl0~ 3 


7T 


±1 


11,12 


.1764959 


0.0 


1.292250xl0~ 2 


7T 


±1 


13 


.1728539 


0.0 


1.158930xl0~ 2 


7T 






TABLE I. Eigenvalues below TV (A). The first 13 eigenval- 
ues Aj and amplitudes Aj for A = — .5,T = .10 and N = 12. 
The value of A is 26.67983. 



As we increase the temperature the eigenvalues remain 
real and the amplitudes all remain negative until we reach 
a temperature TV (A) where a most dramatic effect oc- 
curs. Namely vast numbers of the amplitudes all vanish 
and the eigenvalues with vanishing amplitudes occur in 
degenerate sets in which all allowed values of k occur for 
N = 2(mod4). (For N = 0(mod4) the eigenvalue with 
k = N/A is degenerate but the amplitude does not van- 
ish). As N — > oo there are an infinite number of these 
degenerate sets. For A = —.5 and N = 12 we find that 
TV(-.5) is slightly below .137832224. In table 2 we give 



III. COMPLEX EIGENVALUES 

As the temperature is increased above the universal 
vanishing temperature TV (A) we very soon reach a tem- 
perature where two real eigenvalues collide and become 
complex conjugate pairs. This happens first with eigen- 
values which are very small. As T increases beyond the 
temperature of first collision Tc(A) successively larger 
eigenvalues collide. At each collision the amplitudes di- 
verge however the correlation functions remain finite due 
to a perfect phase difference of it at the point of level col- 
lision. At the temperature of collision the relevant expo- 
nential decay terms in (||) are replaced by ^4n(Aj/Ao)' i_1 . 

We illustrate the collision of eigenvalues in table 3 
where for A = —.5 and N = 12 we give the first 13 
eigenvalues and amplitudes at T = .151493876679 which 
is slightly above the collision temperature for the pairs 
of eigenvalues 4,5 and 6,7. Because the temperature 
is above collision the eigenvalues and amplitudes of the 
eigenvalues 4,5,6,7 are complex and all four combinations 
of the signs of these two phases are to be used. It is clear 
from the table that the amplitudes are diverging and the 
change in phase of the amplitude from real to 7r/2 sug- 
gests that this divergence must be a square root. Further 
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analysis of the behavior of the amplitudes in the vicin- 
ity of the collision shows that this is indeed the correct 
behavior. 



V. THE COMPLETION OF THE CROSSOVER AT 
THE UPPER CROSSOVER TEMPERATURE 

Tu(A) 



j 


IVAo| 


\j phase 


1^1 


Aj phase 


k 


1,2 


.4481334 


0.0 


4.941102x10^ 


7T 


±1 


3 


.2918946 


7T 


7.509636xl0~ 3 


7T 





4-7 


.1909054 


±1.1899 x 10~ 6 


4.245xl0+ 3 


±1.5708 


±2 


8.9 


.1746361 


0.0 


5.837584xl0~ 3 


0.0 


±1 


10 


.1741270 


0.0 


4.018312xl0~ 3 


0.0 





11 


.1561094 


0.0 


5.198129xl0~ 2 


7T 


3 


12,13 


.1180399 


7T 


3.785983xl0~ 3 


7T 


±1 



TABLE III. A collision of eigenvalues below Tl(A). 
The first 13 eigenvalues Aj and amplitudes Aj for 
A = -.5,T = 0.151493876679 and N = 12. The value of 
Ao is 8.224613. The order of eigenvalues is the same as in ta- 
ble 2. The phases and the magnitude of \Aj\ for eigenvalues 
4-7 have only been determined to four places because of the 
numerical instabilities resulting from the divergence of Aj. 



IV. THE ONSET OF OSCILLATIONS AT THE 
LOWER CROSSOVER TEMPERATURE T L (A) 



In tables 2,3 and 4 we note that there is a positive 
eigenvalue with k = which has a positive amplitude. 
At the lower crossover temperature this eigenvalue does 
not have the largest magnitude and has j = 6 in table 4 
lying below the complex conjugate pairs 1,2 and 3,4 and 
the negative eigenvalue with the negative amplitude at 
j=5. However, as T is increased from Tl(A) the mag- 
nitude of this eigenvalue increases relative to the magni- 
tude of the leading complex conjugate eigenvalue and at 
a temperature Ty(A) the magnitudes of these eigenval- 
ues cross. For A = — .5 and N = 12 this occurs very near 
T=.3611. In table 5 we give the first 13 eigenvalues at 
a slightly higher temperature where the real eigenvalue 
has the largest magnitude. Such a crossing without col- 
lision is possible because the eigenvalues lie in different 
sectors k. We call the temperature where this crossing oc- 
curs the upper crossover temperature. The dependence 
of Tu{— .5) on TV is given in table 6. This N dependence 
is well fit by Tu(-.5;N) = .4328 - .86 /N and this has 
been used to extrapolate to N — > oo. 



As the temperature is increased further we reach a 
temperature TV (A) where the two eigenvalues with the 
largest magnitudes collide and produce a complex conju- 
gate pair. These eigenvalues are in the sectors k = ±1. 
For A = — .5 and N = 12 this occurs slightly below 
T = .22782901. In table 4 we give the largest 13 eigen- 
values at this temperature. The dependence of Tl(— .5) 



j 


IVAol 


Aj phase 




Aj phase 


k 


1-4 


.2388172 


±1.155xl0 -4 


1.9078 xl0+ 2 


±1.570925 


±1 


5 


.1815640 


7T 


7.424416xl0~ 3 


7T 





6 


.1779166 


0.0 


2.511491xl0~ 2 


0.0 





7-10 


.1347183 


±.4576581 


2.968286xl0~ 2 


±2.561349 


±2 


11 


.0903745 


0.0 


4.730570 xl0~ 2 


7T 


3 


12,13 


.0641918 


7T 


2.135921xl0~ 3 


7T 


±1 



TABLE IV. Eigenvalues at T L (A). The first 13 eigenvalues 
Aj and amplitudes Aj for A = -.5,T = 0.22782901. and 
N = 12. The value of A is 4.173698. 



on N is given in table 6. This N dependence is well fit 
by T L (-.5;N) = .2015 ± 3.76/iV 2 and this has been used 
to extrapolate to N — > oo. 

At the lower crossover temperature the leading asymp- 
totic behavior of the correlation is 

S z (n;T L (A), A) ~ ^(A/Aq)"" 1 . (6) 

As the temperature increases further the leading behavior 
of the correlation is oscillatory with a wavelength which 
decreases as T increases. 



J 


34 


|A,/Ao| 


Aj phase 


\A,\ 


Aj phase 


k 


1 


6 


.1862180 


0.0 


6.113537e-02 


0.0 





2-5 


1-4 


.1858948 


±.5968364 


5.134231e-02 


±2.118513 


±1 


6-9 


7-10 


.0863447 


±.7691693 


2.144874e-02 


±2.723044 


±2 


10 


5 


.0850952 


7T 


3.715138e-03 


7T 





11 


11 


.0428430 


0.0 


3.407129e-02 


7T 


3 


12-13 


12-13 


.0275197 


7T 


6.337585e-04 


7T 


±1 



TABLE V. Eigenvalues at T V (A). The first 13 eigenvalues 
Aj and amplitudes Aj for A = —.5, T = .3620 and N = 12. 
The value of Ao is 2.679573. The column j4 is the value of j 
of the corresponding eigenvalue in table 4. 



For all temperatures above the upper crossover tem- 
perature the leading eigenvalue is positive and has a pos- 
itive amplitude. Consequently for large n the correlation 
will be monotonic and positive. 



N 


T v 


T L 


T v 


10 


0.1653900 


0.2391 


0.3463 


12 


0.1378322 


0.2278 


0.3611 


14 


0.1181419 


0.2207 


0.3716 


16 


0.1033742 


0.2162 


0.3791 


oo 


0.0 


0.2015 


0.4328 



TABLE VI. The numerically determined values of the tem- 
peratures T V (A),T L (A) and T V (A) with A = -.5 for various 
values of N. 
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VI. CONCLUSIONS 

In [ffl we found that there is a range of temperatures in 
which the correlation function S z (n;T, A) changes sign 
from negative to positive as n is increased and in we 
discussed level crossing as a mechanism for this effect. 
The studies of the eigenvalues of the quantum trans- 
fer matrix presented here demonstrate that instead of 
the single zero seen in the previous finite size studies 
there is a range of temperatures from Tl(A) to Tjj(A) 
where the correlation function has an oscillatory behav- 
ior with a temperature dependent wavelength. This phe- 
nomenon has not been seen before in any study of the 
XXZ chain and can only occur because the non hermi- 
tian quantum transfer matrix can have complex eigen- 
values. These eigenvalues can be studied analytically 
by using the Bethe-like equations derived in How- 
ever, the necessity of considering complex eigenvalues re- 
quires considerations which have not been encountered 
in previous studies of related Bethe Ansatz equations. 
As stated previously Q we believe that the quantum- 
classical crossover which is seen in the quantum trans- 
fer matrix is deeply connected with the bound states at 
T = seen in the spectrum of the Hamiltonian (|J). An- 
alytic investigations of these questions will be pursued 
elsewhere. 
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